NPS  ARCHIVE 
1967 
MILLER,  A. 


PROPAGATION  OF  'TWO  SIMPLE  ACOUSTIC 
TRANSIENTS  IN  AH  l&OVELGCITY  LAYER  WIT 
PSOTCTLY.REFLECTING   BOUNDARIES 

ALAN  CHARLES  MILLER 
&  GORGE  FREDERICK  WENDT 


Hull 


mi 


WWISi— 

WmwmMMM 


tJm* «« 

iJURranuiwI 


ilSi  '''■''"^';V'  '■■■' 

I  I  '^' ll  '■  '■ ;- '■"■■■  i ■■'■  \l ; ■'■■'. ■'- ■  ■  ■  "    : '  ■  ' fi 


PROPAGATION  OF  TWO  SIMPLE  ACOUSTIC  TRANSIENTS  IN  AN 

ISOVELOCITY  LAYER  WITH  PERFECTLY -REFLECTING  BOUNDARIES 

by 

Alan  Charles  Miller 
Lieutenant,  United  States  Navy 
B.S.  Chemical  Engineering,  Iowa  State  College,  1959 

and 

George  Frederick  Wendt 
Lieutenant,  United  States  Navy 
B.S.C.E.,  Ohio  University,  1959 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1967 


ABSTRACT 
The   propagation  of    simple   acoustic   transients   in   an 
isovelocity   water   layer   'was   investigated  for   the   purpose   of 
studying  a   relationship  between   the   propagation  of   transients 
in   ducts   and   in  layers.      The   Laplace   transform    method   was 
used   to   obtain  approximate   solutions   for   the   acoustic 
pressure,    particle  velocity,    and  particle   displacement  result- 
ing from   a   step-function  input   in  velocity  and   the   particle 
velocity  resulting  from   a  gated   sine-wave   input   in  velocity. 
Computer   programs   were   written  to    evaluate   and  graph 
the   resulting  waveforms.      The   waveforms   resulting  from 
step-function  and   gated   sine-wave   inputs   in  velocity   were 
observed  and   compared   with   the   predicted   waveforms.      With 
the  use   of   Mylar    transducers,    good   correlation  was   obtained 
in  the   region   of   validity  of    the   theoretical   solutions.      It 
was   found   that   previously  published   solutions   for   the 
propagation  of    transients   in  ducts  became   the   leading  terms 
of    the   associated   solutions   for   the   layer. 
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1.  INTRODUCTION 

While    most   investigations    of    the   propagation   of    energy 
in   waveguides   have   been   done   for    electromagnetic    cases, 
some   have    considered   the   propagation   of   pulses    of   acoustical 
energy   in   waveguides^     '  2>3  tk)  ,       The   dispersion   of    an   ultra- 
sonic  pulse   propagated   through   a   waveguide   of    circular    cross- 
section   filled   with    water   has   been   used   for    the   approximate 

determination   of    the   frequency   spectrum   of    the   original 

(  1) 

pulse         .       It  has   been   showi    that    the    dispersion   of    a 

rectangular   pulse   along   an   acoustic    waveguide   of    rectangular 
cross-section   can  be   used   to    calibrate   the   waveguide   so    it 
can  be   used   to    predict   the    effects    of    similar    dispersion   on 

the   transmission   of    pulses    of    other    shapes    that   are   not 

(  2) 

amenable   to    actual    computation     '  .       Pulse-compression 

studies    in  an   acoustic   waveguide    show    that   a   frequency 
modulated   input   pulse    can  be    compressed   into   an   output   pulse 
of    shorter    duration   and   of   higher   peak-amplitude^    .       In 
addition    to    the    waveguide    experiments,    investigations    of 
transients    in   layered    media  have  been    made      '       .      W.    C. 
Knudsen^ J/    investigated    the   propagation   of    a   pressure 
transient   in  a    surface   layer    of    liquid    which   lay   over   a   deep 
bottom    layer    of    liquid.       Many   of    these    studies   have   used 
approximation   techniques    such   as    the    method   of    stationary 
phase.      A   summary   of    these    techniques    useful   for    calculating 
the   resultant   propagated   waveform   has   been   presented  by 
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«J .    R.    Wait^7).       The   original   work   of    Rayleigh   on   scatter- 
ing  and   radiation   of   harmonic    disturbances   by   small 
cylinders    was    extended   to    transient   disturbances   by   «J .    W. 

I  a) 

Miles    using   the   Laplace    transform^0'.       R.    B.    Barakat 
used   the  Laplace    transform    to   investigate   the   radiation   of 
acoustic   pulses   of   various    shapes   from   an  infinitely   long 
circular    cylinderV7/. 

Recently,    exact   solutions   for   the   propagation   of    simple 
transients   in  waveguides    (step-function  and   gated   sine-wave 

(  )      1  n     1  -I  \ 

inputs)    have  been  published      '       '         .      Some   work  has   also 
been   done   which   suggests   a   correlation  between   these   exact 

solutions   in   ducts   and   solutions   for   the   propagation  of 

(  12) 

acoustic    transients    in  isovelocity   layers 

The   purpose   of    this   research   is    to   investigate   the 
region   of    validity   of    the    suggested    correlation   and   to    see 
what,    if   any,    restrictions    must  be   placed  on   the   resulting 
solutions    describing   the   transient   response   of    the   layer.       To 

facilitate    this    comparison,    those   results    of    the   propagation 

{ )     1  P  ) 
of    transients   in   ducts      '  which   will  be   of   value   are 

summarized   in  Appendix  A. 
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2.  THEORY 

Development   of    Solutions 
The   response    of    an   isovelocity   layer    to    simple   transient 
excitations    will  be   developed   using   the   following   assumptions: 

(a)  Cylindrical   velocity-source    of    radius    a. 

(b)  Constant   layer   depth,    d. 

(c)  Perfectly-Reflecting   upper    and   lower   boundaries. 


An   Isovelocity  Layer 
Figure   2. 1 
For   radially-symmetric   propagation   in   the   layer   indicated 
in   Figure    2.1,    the   wave    equation   for    the   acoustic   velocity 
potential    C£    can  be    written   in    cylindrical    coordinates    as 


(&+H.*|ji-i£)$<v#-° 


(2.1) 


where   r   is    the   radial   direction  and    z    is    the    transverse 


direction. 
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The    separation   of   variables 

leads    to    two   partial   differencial   equations:    one   for   a 
"transverse"    function    £^  which   depends    only   on   the   height 
z    of    the   layer 


[$-par:K<»-o, 


(2.3) 


and   one  for   a   "radial"    portion  0n   of    the   acoustic   velocity 
potential, 


C6>t* -(#-*&><« -o 


(2.4) 


which   depends    on  both   time,    t,    and   the   distance  from    the 
origin,    r.       The   value   of    the    separation   constant  COn  depends 
on   the  nature   of    the  boundaries   at    z   =   0    and   d.       For 
example,    if   we  have   a   pressure-release   top   and  bottom, 
then  p   =  -P  2i     =0    at    z   =    0    and   d  and   Eq.    2.3    has    the 
solution 

c 

where  o>n   =  n~ITC/j      •       Furthermore,    since   the    source   is 
assumed    symmetric   with   respect   to    the   plane    z   =    d/2.just 
those    modes   having  pressure   antinodes    at   d/2    will  be    excited , 
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Thus    only   odd   values    of   n   are   allowed.       Therefore    Eq. 
2.5    can  be   rewritten   as 


(2.6) 


Since    this    transverse   function   affects    the   form    of    Eq.    2.4 
only   through   the   allowed   values    of  CO    ,    it   is    convenient   to 
establish   a   set   of    reduced   acoustic   variables    which    suppresses 
the   function    ^^^3) 

The    well-known   acoustic   relations 


(2.7b) 


and 


^rn!      |rI„  (2-7=) 


where   n    is    the    mode   number,    p      the   acoustic   pressure   for 
the   nth    mode,    u      the   particle   velocity   for    the   nth    mode, 

and   u        the   radial   particle   velocity   for    the   nth    mode    can  be 

rn  *  J 

used   with    Eq.    2.2    to    show    that 


■T„»-c  *»<})&4,tr,t).  (2-8) 


bt 
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The    "fragment    "    p  of    the   acoustic   pressure   for   the   nth 

mode   is    defined   as 


•f,„=  %»  =  -f# 


where   0      =   0      (r,t)    is    the   radial   portion   of    the   velocity- 
potential   for   the   nth    mode   which   -was   previously   defined  by 
Eq.    2.2.       Combination   of    Eqs.2.2   and   2.7c   allows    the 
definition   of    the    "fragment"    of    the   radial   particle-velocity 
for    the   nth    mode, 


^-5     ""/*,«    S    &*•*  (210) 


Since    the    transverse   function    Z      (z)    is   independent   of 
time,    the   radial   particle-displacement  for   the   nth    mode    may 
be   obtained  by   integration   of    u        with   respect   to    time.       This 
enables    the    "fragment"    of    the   radial   particle-displacement  for 
the   nth    mode   to   be   written  as 

These   reduced   acoustic   variables   or    "fragments"    0n,      f\p        , 

p      ,    and     ^  will  be   used   throughout   the   remainder   of    this 

rn  0    rn 

development. 
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Since    we   are    considering   the   propagation   of    transients, 
Eq.    2.1+t    a   differential    equation   in   two    unknowns,    cannot  be 
solved  by   separation   of    variables.       The   Laplace    transform 
method   will   be    used    to   obtain   the   desired   solutions. 

The   Laplace    transform    of    Eq.    2.1+    with   quiescent 
initial    conditions    results    in   a   differential    equation   in   one 
unknown   which    may   be   written  as 


[&*it -(&*)!*. -° 


(2.12) 


where    the   bar    over    the   function   denotes    the   Laplace   trans- 
form   of    that   function: 


The    change    of    variable 


-   .rV^+ 


■» 


>C    =•     rV^+^n  (2.13) 


casts    Eq.    2.12    into    the   form 


(2.14) 

The    solutions    to    Eq.    2.8   are   I       (ytc  )    and    K      (  Xc )  ,    the 
modified   Bessel   functions    of    the   first   and   third   kind*    ^    . 
Conservation   of    energy   requires     0^     to    go    to    zero    as    s-*oO, 
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so   that    the   only   acceptable   solution   of    this    differential 
equation   is    K      {/JL. )  •       The   Laplace    transform    of   0n    can 
now   be    expressed   as 


£    =   CKjjA 


(2.15) 


where   C   is   an  amplitude    constant  to   be   determined  from 
the   boundary    conditions. 

With    the   use   of    the   relationship   expressed  by   Eq.    2.7c, 

the  Laplace   transform    of  AT         becomes 

rn 


Kn  -i^cy^te*  c«r™      (2.i6) 


The   known  boundary   condition  at   the   source,    located   at 
r   =   a,    for   propagation  in   the   nth   normal   mode, 


(2.17) 


has    the    transform 


/\T     K1)    =  Y     .  (2.18) 
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Equations   2.16    and   2.18    allow    evaluation   of    C, 


c    - 


so    that     Of"  can  be    expressed   as 

y     rn 


(2.19) 


(2.20) 


The    substitution   for    C   in   Eq.    2.15    allows  ^n   to   be   written 


as 


The    expression  for    the   Laplace   transform    of    the   fragment 
of    acoustic    pressure,    p      ,    can  be   obtained   from    Eq.    2.21 
if    we   recall    that   Jf      =    —P  *+JlT    corresponds    to    -jo      s-^S/i, 

in   the   transform    domain.       Thus, 


js   _    rest  K.(glftw') 

urn  v  "— * 


(2.22) 
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The   Laplace   transform    of    the   particle-displacement 
fragment,   \  ,    can  be   obtained  from    /V"rn  by   either 

integration   "with   respect   to    t   in   the   time   domain   or    multi- 
plication by    1/s    in   the    transform    domain.       The   result   is: 

Inverse   transforms   of    Eq.    2.20    -   2.23   have   not  been 
tabulated,    so   approximation   methods    must  be   used. 

If  yL^  (Eq.    2.13)    is   greater   than  unity  the   asymptotic 


form 


K»(^=i/i5  «-"A('+  ifir*— ) 


(2.24) 


( 12    13 ) 
may  be   used         '         .       The   requirement  M  >  1    can  be 

satisfied  by  restricting   c/au)      to   values   less   than  unity. 

Substitution   of    Eq.    2.2i+   into    Eq.    2.20    yields 


tfrn=. "\Uj7cfeggg    g*pC-efsWJ(l+fl£V^a?)       (2.25) 
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which    simplifies    to 


(2.26) 


where   T    =   r-a   is    the    time-of -flight   of    the    signal   from 
c 

source    to   receiver.       Since    c/acO  n   is   restricted   to    values 
less    than   unity,    the   denominator    can  be    expressed   as 


(l  +  gfeW?) 


so    that   we   have 


l~    Sf^v^T        5  (2-27) 


v^-VjIF  «*L-t«W]  t+  «I^55?]  t~  fife] 


(2.28) 


which    simplifies    to 


1?  ^m^  ^P[-T\f7H3?]  -j  3c  (j.^  CXf{-T^+^]  (  2  .  29 ) 


8aAkV^rF 


2 
where  terms  of  order  (  c/a  to  )   have  been  dropped. 

n 


23 


With    similar    manipulations    the   other    desired   acoustic 
variables    can  be   obtained   from    Eqs.    2.21   -   2.23: 

(?0n=  -Kg   ^pEW^?J  +L  3c1-  /i+A\  <2xp£t£^?1  (2.30) 


(2.3D 


8o_\|s^ 


Equations    2.29   -   2.32   approximate   the   Laplace   transforms 

of    the   desired   acoustic   variables   for    the   boundary   condition 

Va(t)    at   r    =   a   with   the   restriction    c/a  CO        ^.     1. 

n 

Comparison   of    Eq.    2.29    "with    Eqs.    A.  7    and  A. 9    shows 
that   the    terms    of    Eq.    2.29   have   known   inversions.       For   a 
step   input,    for   which  V      =      /s    the   inverse   transform   for 
the   fragment   of    the   radial   particle-velocity    can  be   written 
down   directly  from    the   known   results    of    the   duct    (see 
Appendix  A. 

(2.33) 
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J 


where  ARG   =  OJ      \/    t   -T      and   it   is   always    understood   that 


all    solutions    vanish   for    t  £  T . 

The    inverse   transforms    of    the    other   acoustic   quantities 
can  be   found    similarly   by   using   the   known   transform   pairs 
from    the   duct    case   and    the    additional   transform   pairs 
developed   in  Appendix  A.       The   results    are 


r  =6 


\Ja_  rV>~ 


(A* 


.lm-» 


06 


£(£P'Tm  -/H1^  £(&HJ>* 


(2.34) 


|£f«,»r 


(2.35) 


nM!'-Wl-f>-s.^imfj^ 


(2.36) 


The    quantity   t'    =    t   -   T    appearing   in    Eqs .    2.33    -    2.36    and 

2  2 

in  ARG   as(t-T)(t   +    T)=t      -    T      represents    the    time 

which   has    elapsed   since    the    signal   traveling   with   velocity    c 

first   arrived   at   the   position   r.       It    will   be   referred   to    as 

the   delayed    time. 

Equations    2.33    -    2.36    are    the    solutions   for    the    transient 
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response   of    an  isovelocity   layer    to    a   step   input   in   velocity 
at  r   =   a   and   t   =   0 . 

Similar   procedures    may  be   used   to    obtain   the  fragment 
of    the   radial   particle-velocity  for   a   gated    sine-wave   input 
in   velocity   for    which     Y^-  "^7\%   * 


^^^-(j/^rtRfcO     (2.37) 
_  — ' 


where  l/i(fUQ)rt£6]    =   Lommel   function   of    two   arguments      ^' 


VHLk^ 


and        AXTi  =-    <o(-£-TY  /  *  <%/c)# 


Other    acoustic    quantities    could   not   be    completely   solved   in 
this    case.       Details    are    shown    in   Section  A-3    of   Appendix  A, 

Comparison   of    the   above    solutions   with    the    solutions 
for    the    duct   listed   in  Appendix  A   shows    that   in   all    cases 
the    solutions   previously   obtained   for    the   duct   become    the 
leading   terms    in   the    solutions    developed  for    the   layer. 

A   physical   understanding   of    this    can  be   gained  by 
comparing   propagation   processes   in  a   duct   with   those   in   a 
layer.      While   a   duct   presents    a    constant    cross-section   to 
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the   propagating   wave   form   for   all    time,    a   transient 
propagating   outward   in   a   layer   is   presented   with   an 
effective    cross-sectional   area   which   increases    as    r.       How- 
ever,   for    sufficiently   large   r    the    cross-sectional   area   will 
not    change   appreciably   as    the    waveform    propagates    through 
unit   distance.       In   the   limit,    therefore,    we    expect    the 
radial   divergence    contained  by   the   term  — •  =J-  in    Eq.    2.4 

to   become   insignificant.       In   this   limit,    propagation   processes 
can  be    described  by   an   equation   virtually   the    same    as    Eq.    A3 
which   describes    propagation   in   a   duct. 

For    this    reason   it   is    not    surprising   that   the   solutions 
for    the   duct   become    the    leading   terms    of    the   layer    solutions. 
The    remaining   terms    in   the   layer    solutions    are    correction 
terms    necessitated  by   the   fact   that   the   layer    does    not 
actually   present   a    constant    cross-section   to    the   propagating 
wave   form. 

The    solutions   for    the   transient   response    of    an 
isovelocity   layer   have   been   developed   with   the   restriction 
c/a  u>       <.  1.       Since   the    cutoff   frequency   is   given  by 
(jj    =nTrc-/e/        »    this    restriction   places    a   limit   on   the 
relation  between   the   depth   of    the   layer    and   the    size   of    the 
source,    d/a  y  k>TT    .       If    it   is    desired   to   remove    this   restric- 
tion  a   different    expansion   of    the  K,.(/x_\     's    involved    must  be 
used. 
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Investigation   of    Solutions 
The    solutions    developed   for    the   layer   have  been  obtained 
for   a    single    mode.       In   determining   the    contributions    to   the 
solutions   from    other    modes,    attention   must  be   given   to 
the    transverse  behavior    of    the    transmitted   and   received 
signals.       The  boundary    condition   on  velocity   at   the    source 
at   r   =   a   says   that  V    (t)    vanishes   for    z    <0    and     z      >d  and 
is    constant   for    0  <,    z  <^  d.       This    condition   can  be    met   by 
representing  V    (t)    by   a    spatial   Fourier-superposition   of 

[u(z)    -   u(z-d)J    in   terms   of    the   normal    modes    across    the 
face   of    the    source: 

1 

Notice    that   the   higher    modes    contribute   as    ^     «  for    this 

input.       These  higher    modes   have   positive   and   negative 
portions   in   their    waveforms    over    the   depth   of    the   layer. 
These   positive   and   negative   portions    will   partially    cancel 

■      l|«4   '.III* 

over  the  face  of  a  receiver  which  extends  over  the  layer 

depth.   In  effect  the  receiver  will  average  the  received 

(1) 

signal   over   the   area   of   its   face         : 


1/L    » 
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Evaluation   of    the   above   integral   gives 


<ad> 


nr     =.     8d     y  /-L\L  (2.1,0) 

.    i   a 
Therefore    each    mode    contributes   as   (r— -J  or    1,    1/9, 

1/25,    1/49,    for    the   first   four    excited    modes. 

A    computer   program    was    written   to    evaluate   the 
effects    of    these   higher    modes    on   the    solutions    developed. 
Computer   predictions    considering    (a)    the   fundamental    mode 
alone   and    (b)    the   fundamental    mode   plus    the   next   four   non- 
zero  higher    modes,    were   plotted   versus    delayed-time. 
Comparison   of    these    two   predictions    as    presented   graphically 
by   the    computer    showed    changes    in   the    amplitude    of    less 
than    1    per    cent   and    shifts   in   axis    crossings    of    less    than 
0.01    local      wavelengths.       Since    these    effects    are    negligible 
compared    to    the    experimental   errors    encountered   in   this 
research,    it   was    sufficient   to    retain   only   the   fundamental 
mode    ( n   =    1)  . 

Examination   of    Eq.    2.36    shows    that   it    consists    of 
two    terms.       The   first    (main)    term   is   from    the   duct 
solution   and    the    second    (correction)    term    is   from    the   layer 
effect.      An    estimate   of    the    effect    of    the    correction   term 
can  be    made  by    combining   the   leading   terms    of    the   indicated 
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sums.       Equation   2.36    can  be   rewritten   as 


Under    the   assumption   t'    =    t-T  <^    T,    which   limits    considera- 
tion to   the   leading  portion  of    the   transient   waveform,    we 
have   tftT,    and    Eq.    2.41   becomes 


(2.42) 


w 


here     ©        =  ARG  =     04^  t' Z.T' 

and        B»  -*&,(!-*>(&)*. 

Substitution   of    trigonometric   approximations   of    the   Bessel 

( 13) 
functions  in    Eq.    2.42   gives 


^"^s  i  cU(e-^.  w«e-y).i*- 


43) 


By   combining  the    cosine    terms    of    Eq.    2.43>    we  find   the 
expression  for    the  fragment   of   radial  particle   displacement   to 
be 
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where         £f         =   tan"   B(t). 

Equation   2.1+1+   shows    that    the    correction   term   in   the 
layer    solution   introduces    a   phase   angle      Q        that   was   not 
present   in    the    duct    solution.       The    zeroes    of    the   right  hand 
side    of    Eq.    2.1+1+   give    the    times   for    which   the   particle 
displacement   vanishes.       It  has   been   shown   that    the   duct 
solution  becomes    the   leading   term    of    the   layer    solution; 
therefore,    the    zeroes    of    the   first    term    of    Eq.    2.1+3    give    the 
nulls,    or   axis    crossings,    of    the   particle   displacement   in  a 
duct.       If    these    times    are    evaluated,    the    shift   in   time    of 
associated   axis    crossings   between   the   layer   and    the   duct 
solutions    can  be   predicted. 

The    zeroes    of    the    layer    solution,    designated  by     £*  , 

are   obtained    when 


COs(6-$J?+)S)   =  0;  (2.45) 


or 


9-^^  --     ^d)Tr3m.i)J)Vtc.      (2.Mi) 
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The    zeroes    of    the   duct    solution   will  be    designated  by     7~ 
They   are   obtained   when 

(2.47) 


cos(e>-\jr\  &  o} 


or    when 


(^-^")  -    (€£22z'V     ms^2#BJe^e.       (2.48) 


Since   we  have   t'  <  <    T,    c/a  to       <,     1,    and  r  r^r  a,    so    that 


-W  V  -  *tt  =  1^(1-^^  <^  1  ^      (8iW) 


a   small-angle   approximation   can  be   applied   with    the   result 


+o.n  X  =  K^t1  (2-50) 


where 


1  u.  o.^  V.1     v-  ) 


\|3T 
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Substituation   of    this    expression  for      Q  into    Eq.    2.1+6    and 

evaluation   of      "       in   Eqs.    2 .  i+6    and   2.1+8    gives 


and 


("v^V^  -(2^i^)7r. 


Manipulation   of    Eqs.    2.51    and   2.52    reveals    that   the   ratio 


-      / 

#vy  ^m       can   be   expressed   as 


r~  1 


»v\ 


Thus,   ■— V  '0-<»V»)  is    a    convenient   parameter    to    use   in 

predicting  the  shift   in   axis    crossings. 

This   result   is   approximate,    since    only   the   first   terms 
of    the    sums   in    Eq.    2.36    were   used   to    develop   it.     To 
evaluate    the   accuracy   of    this   approximate   result,    computer 
programs    were   designed   to    evaluate   the   predicted   waveforms 
of    Eqs.    2.33    and   2.35   -    2.37.       They   were   written   so    that 
the   main   terms    and   the    correction   terms    of    these    equations 
could  be    evaluated    separately.       This   permitted    comparison   of 
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the   layer    solutions    and   the   duct   solutions    and   showed   the 

behavior    of    the    correction   term   as    a  function   of    c/a  Cd 

n 

These   programs,    one   for    the    step   input    (Appendix  A)    and 
one   for    the   gated    sine-wave   input    (Appendix   D)    consider 
up   to    11    terms    of    the    sums    involved.       The   programs   give 
printed   and   graphical   output    of    the    main   terms,    the 
correction   terms,    and   the    complete    expressions    of    the 
equations    mentioned   above.       Figure    2.2    shows    the   velocity- 
fragment   predicted  by   the    computer.       Figures    2.3    -   2.5 
show    the    complete    expression,    the    main   term,    and   the 
correction   term    of    the    displacement   fragment   for    several 
values    of    c/a  LJ      . 

Information   about   the    envelope   of    the   displacement 
fragment   waveform    can   also   be   obtained   from    Eq.    2.i+4 


which    gives    the   amplitude   as         \J\"Tl3("t)  .       Under    the 

assumption   t'    <<     T,    B(t)    becomes 


Bto 


(2.54) 


and 


^■■k&mTt 
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Therefore,    the   amplitude   of    the   initial   portion   of    the 
waveform    under    consideration   is 

The   approximate    solutions   for    the    propagation   of    simple 
transients   in   an   isovelocity   layer   have   been   obtained,    and 
computer   programs    to    evaluate    these    solutions    and   predict 
the   axis-crossings    and   amplitude-modulation   in   the   layer 
have   been   written. 
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(a)     Step   input  with  c/ac^s  O.69,    rs6l.5  cm, 
a  r  2.54  cm,   and  dr^.^  cm. 
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(b)     Gatad  «ina  wav«  input  with  c/aH,Bl»00 

r«6l#5  cm,  am  2.54  cm,   d*7.97  cm,  cj  -  585  IcHz. 


Fragment  of  radial  particl*  velocity,  v^n  ,  v«.  dalayad  tim«,  t». 

Figur*  2.2 
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(a)  Complete;  sxpression  (layer  solution) 
Main  term  (duct  solution) —  — 
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(b)  Layer  correction  term 


Fragment  of  the  radiel  particle  displacement,^**  vs*  delayed  time, 
t1,  as  predicted  by  the-  computer  program.  c/auH-0.35,  r*  6l»5  cm, 
a  s2.54  cm,  d-  2.80  cm. 


Figure  2.5 
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(a)  Complete  expression  (layer  solution) 

Main  term  (duct  solution) 


(b)  Layer  correction  term 


Fragment  or  the  radial  particle  displacement,  So»>  vo*  delayed  time, 
t1,  as  predicted  by  the  computer  program.  c/aon=1.00,  r-  6I.5  cm, 
a  =  2.54  cm,  d^7»97  cm. 


Figure  2.4 
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Complete  o;:prc  oion  (layer  solution) 


Main  term  (duct  solution) —  -  - 


Fra0ment  of  the  radial  particlo  displacement,  3m  >  V8«  delayed  time, 


t',  as  predicted  by  the  computer  program, 
a-  1.00  cm,  d  =  9.^2  cm. 

Figure  2.J?a 


c/au  =  5'QO,  r«6l»5  cm, 
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Layer  correction  term 


Fragment  of  the  radial  particle  displacement,^,  vs.  delayed  time, 
t',  as  predicted  by  the  computer  program.  c/au>H=3»00»  r-61.5  cm, 
a=  1.00  cm,  d  =  9.42  cm. 

Figure  2. 5b 
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3.  EXPERIMENTAL   APPARATUS 

The    experimental   apparatus    used   in   this    research 
consisted   of    a   tank,    several   transmitters   and   receivers, 
and   the    associated   electronics    equipment.       The   laboratory 
setup   is    shown   in   Fig.    3.1. 
Tank 

The   water    layer    was    contained   in   a   tank    constructed 
of   3/4   inch    marine   plywood   which    was   30    cm   high,    1    m    wide, 
and    1.5    m    long.       The   plywood   was    coated    with   Varathene 
to    provide   waterproofing,    and    lined   with    1/2+   inch 
polyethylene   foam.       This    provided   the   desired   pressure- 
release  bottom   boundary   for    the   layer.       The    water-air 
interface   provided    the    top   pressure   release    surface.      With 
the    sides    and   ends    of    the    tank   also    lined   with   polyethylene, 
the   tank    could  be    excited   as    a   pressure-release    cavity   which 
served   to    accurately   determine    the    depth   of    the   layer. 
This    technique   is    described   in    more   detail   in   section   i+. 

Cylindrical    sources    could   not  be   used  because    the   tank 
was    too    small.       Transmitters    were    therefore   built   in   the 
shape   of    semi-cylinders    and    mounted   on   a   rigid   plate.      A 
neoprene-coated    steel   plate,    3/8   inch    thick,    was    made    to 
fit   in   one    end   of    the    tank   to    act   as    the   rigid  boundary. 
This   rigid    surface    served   as    an   acoustical    mirror,    so    that 
the   reflected   pressure   appeared   to   be    emitted  from   an 
image    source    (an   identical    semi-cylinder   back-to-back   with 

kl 
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the   real    source).       Under    this    condition    cylindrical    symmetry- 
was    approximated.       Figure   3-2    illustrates    the   geometry 
of    the   layer. 
Transducers 

Mylar    dielectric   transducers    were    selected   to   be   used 
in   the   layer    for    the   following   reasons; 

(a)  Their   broad-band    transmitting   and   receiving 
characteristics    allowed   the    transmission   and   reception   of 
the    transient   excitations    without   introducing   significent 
frequency-sensitive    modulation   of    the    signals. 

(b)  They   could  be    conveniently    constructed   to    the 
required   semi-cylindrical    shapes    and   sizes. 

(c)  Previous    work^       9     ''    has    indicated   that,    when 
driven   with   a   low-impedance   voltage    source,    they   generate 
an   acoustic    signal   for    which   the   normal   particle   velocity   at 
the    source   is   proportional   to    the   input   voltage. 

Two    transmitters    were    constructed   to    the   required 

(  17) 
semi-cylindrical    shape,    using   the   techniques    described   in  , 

one   with   radius    1.0    cm   and   one   with   radius    2.5    cm.       Both 

transmitters    were    12    cm   high    which    corresponded   to    the 

maximum    layer    depth    to   be    encountered.       These   dimensions 

were    chosen   in   order    to    cover    the    desired   range    of   values 

for    c/a  GJ    n. 

A.   Mylar    dielectric   receiver,    identical   in   nature   to    the 

transmitters,    but    with    a   plane,    rectangular   face    measuring 
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3    cm    x   14    cm,    was   used   to    measure    the   averaged    (over 
depth)    radical   particle   velocity   and   displacement   in   the 
layer.       By   proper   adjustment   of   its    electrical   load,    the 
receiver    could   be    made    sensitive   to    either   particle   velocity 
or   displacement.      When   terminated  by   the    1    megohm 
impedance   of    the   voltage   amplifier,    it   was    displacement 
sensitive;    when   shunted   with   an  impedance   of    100    ohms,    it 
became    sensitive    to    velocity.       Transmitters   and   receivers 
are    shown   in   Figure   3«3« 

Two    pressure-sensitive   receivers    were    constructed, 
but   neither   proved   to   be    satisfactory.      A   probe   using   a 
1/8   inch   barium    titanate    cylinder   as    the    sensing   element 
was   built,    but   its    low    sensitivity   resulted   in   an   unacceptable 
signal   to   noise   ratio.      A   2.5    cm    x   2.5    cm    sensing   element 
cut   from    a   0.5    cm    Glennite   disc    was    also    constructed. 
However,    the   resonant  frequency  of    the   element   was 
sufficiently   close    to    the   frequencies    at    the   front    of    the 
received   waveform    so    that   ringing   of    the    element    caused 
severe   distortion   of    the   received   signal.       For    this    reason, 
and  because   it    was   felt    that   the   velocity   and   displacement 
measurements    would   adequately   describe   the   waveforms, 
pressure    measurements    were   not   taken. 
Electronics    Equipment 

Figure   3*4    shows    the   block   diagram    of    the   equipment; 
Table   3*1   identifies    the   items    in   the   figure.       The   oscillator 

kk- 


On  the  left,  a  2.5  cm  radius  backplate  and  finished 
transmitter;  center,  partially  completed  and  finished 
receivers;  on  the  right,  bhe  1.   cm  ra  ':.iu  ■■  transmitter, 

Oransmitters  an  !   sceivers 

E  i  ;u  re   o 
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Table   3-1 

Electronics    Equipment 

OSC  Audio    oscillator,    TS-382    D/U 

TBG  Tone   burst   generator,    General   Radio    type    1396-A 

UPG  Unit   pulse   generator,    General   Radio   type    1217-A 

FREQ  Electronic    counter,    Hewlett   Packard    5233-L 

PWR   AMPPower   amplifier,    Hewlett    Packard   i+67-A 

AMP    1  Voltage   amplifier,    Teletronix   type    1121 

AMP    2  Wide   band    chain   amplifier,    SKL    model   202-D 

RC    FIL      Homemade   variable   low-pass    RC   filter   with 
400    kHz    roll-over. 

SKL    PIL  Electronic  band   pass   filter    SKL    model   302 

SCOPE       Dual   beam    oscilloscope,    Fairchild   model    766-H 

LBN  Load   and  bias   network    (made   for    earlier    work) 

shown   in   Fig.    3-5 


hi 


and   tone   burst   generator    were   used   as    the   initial    signal 
source   for    the   gated    sine-wave   input.       For   the    step   input, 
the   unit   pulse    generator    was   used. 

The   received    signal   was    sent   from    the   receiver 
through    the   load   and   bias   network,    amplifiers,    and   filters 
to    trace    "B"    of    the   dual   beam    oscilloscope.       The   input 
signal    taken  from    the   output   of    the   power   amplifier    was 
displayed   on   trace    "A" .       Measurements    of    the   received 
waveform    axis    crossings   and   amplitudes    were    made   directly 
from    the   oscilloscope   display,    or   photographed   for   later 
analysis.       Triggering   of    the    scope   was   provided  by   gating 
signals   from    either    the    tone  burst   generator   or    the   unit 
pulse    generator.       The   load   and   bias    network   shown   in 
Figure   3-5    served   to    supply   the   300    volt  bias    needed  by 
the   Mylar    transducers    without   allowing   this   voltage    to 
appear    elsewhere   in   the    circuit,    and   to    provide   the   proper 
electrical   load   for    the   receiver    so    that    either   particle 
displacement   or   velocity    could  be   observed.       The   22   ohm 
resistor    in   the    transmitter    circuit   was   needed   to    suppress 
parasitic    oscillations    caused   by   the    transient    electrical 
excitation   of    the    capacitance    of    the    transmitter   and   the 
residual   inductance    of    the   power   amplifier. 

The   rise    time    of    the    system   was    experimentally 
determined   to   be   about   0.5  A\  sec. 
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1+.  EXPERIMENTAL    INVESTIGATION 

The    purpose    of    the    experimental    investigation    "was    to 

(a)    observe   and   analyze    the   received   radial   particle-velocity 

and   displacement    waveforms    in   an   isovelocity   layer;     (b) 

compare    these   waveforms    with   theoretical   waveforms    produced 

as    output   of    computer   programs   based   on   the   theoretical 

developments;    and    (c)    vary   the   parameter    c/a  tJ        over   a 

range   of   values    to    see   if    the    solutions   broke   down    "when 

the   basic   restriction    c/a  o*       <1    was   violated.       The   greater 

n 

part   of    the    experimental   investigation   was    made   using   the 
step-function   input;    the   gated-sine    wave   input   was   used 
only   infrequently.       Experiments    were    conducted   using   the 
procedure   described   in   the   following   paragraphs.       Five  basic 
variables    characterized   each    experiment;    the   input   function 
(gated    sine-'Wave    or    step-function),    the   range   r,    the 
transmitter   radius    a,    the   layer   depth    d,    and   the    speed   of 
sound    c. 
Experimental   Procedure 

Since    c/a  LO        was    a   fundamental   parameter    of    the 
solutions,    the   input   function  and   the   value   of    c/a  ^J       to 
be   used   were   first    chosen  for    each   run.       The   2.5    cm 
transmitter   was    used   for    values    of    c/a  63  n   from    0.35    to 
1.00;    the    1.0    cm    transmitter    was   used  for   values   ranging 

from    1.00    to   3.00.       The   temperature    of    the   water   was 

(  18 ) 

measured   to    determine    the    speed   of    sound    c   in   the   layer v    0/. 
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With   a,    c,    and   the   desired   value   of    c/a  CO        determined 
the   necessary  value   of    LOn     could  be    calculated.      The   tank 
■was    then  filled   to    a    convenient   depth   and   excited   as   a 
pressure-release    cavity   of    known   length   and   "width.       The 
depth   of    water    "was    adjusted   until   the    maximum   resonance 
peak   of    the   fundamental   mode   "was   observed  at   the 
desired  LOn  .       From  CO^  and    the   known   length   and   width   of 
the   tank   the   effective   depth   was    calculated  for   the    computer 
program    which   was  written   with   depth   rather    than    L0n  as 
an   input.       The    end   plate   and   transmitter    were   placed   in 
the   tank,    an   input    signal   was   fed   to    the   transmitter,    and 
the   receiver   was    placed   at   the   range   of    61.5    cm    (used   in 
all    experiments)    by    measuring   the    time-of -flight   of    the 
signal   from    transmitter    to   receiver. 

The   received   waveforms    were   photographed   from    the 
oscilloscope    display.       For    the    step   input,    several   photo- 
graphs  were   taken  of    each   waveform.      An   expanded   time 
scale   was    used   for    photographs    of    the   leading   portion   of 
waveform    where    the   high   frequencies    were    contained;    a 
more    compressed   time    scale   was   used  for    suceeding  photo- 
graphs.      For   a    typical    experiment,    three   photographs    of 
each    waveform    were    made   utilizing,    for    example,    5  /U.sec/ 
cm,    10    yusec/cm,    and    50    yu.sec/cm    time    scales.       For    the 
step   input    experiments,    both   received   particle   displacement 
and   velocity   were   photographed   in    most    cases.       The   photo- 
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graphs    were    then   read   using   dividers   and   a   scale   to    determine 
the   values    of    delayed   time    £*       at   which   axis    crossings 
occured.       The   amplitudes    of    the   peaks    in   the   waveform    were 
measured   in   a    similiar   fashion.       Sample   photographs    are 
shown   in   Fig.    i+.l.       Data   read  from    all   photographs    are 
contained   in  Appendix  B . 

For    each    experimental   run,    a    corresponding   computer 
run   was    made   using   the    same   values   for   r,    a,    c,    and   d 
(and   hence    c/a  to      )    as   inputs    to    the    computer   program. 
As   has   been   shown   in   Section   2,    a    solution  for   particle 
displacement   or    velocity   in   the   layer    can  be   expressed   as 
the    sum    of    the    corresponding   duct   solution   and   a    "layer 
correction."      The    computer   programs    could   therefore   be 
written   so    that   their   output  for    each   run   contained   a 
series    of    graphs    of    the   duct   solution,    layer    solution,    and 
"layer    correction"    waveforms.       The   program   for    the    step 
input   produced   a    series    of    graphs   for   both   particle   dis- 
placement  and   velocity   waveforms.      The  program   for    the 
gated-sine   input /graphed? "the   velacifcyh series   only.       Comparison 
of   graphs    within   a    series   provided   convenient   visual   evidence 
of    the   relative   importance   of    the    "layer    correction"    with 
respect   to   the   duct   solutions    as   a  function   of    the   parameter 
c/a  60  n«      A   portion   of    the    step   input   program    was   used 
to   find   the   axis    crossings    of    the   theoretical   displacement 
waveforms   for    the   duct   and   layer   solutions. 
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Another   portion   of    the   program    evaluated   the    magnitude 
and   delayed   time   of    the   amplitude   peaks    in   the   layer-solution 
displacement   waveform.       In   this    manner    the   theoretical   values 
of    2"         and    L^    ,    and    the    magnitudes    of    times    of    the   peak 
amplitudes    could   then  be    compared   with    the   observed   values 
of  £  rn      and   amplitude   peaks   found   in   the    experiment   with 
the    corresponding   value   of    c/a  CO 
Analysis    of    Data 

Analysis    of    the   response    of    the   layer    to    a    step- 
function   input   proved   to   be    more   involved   than   was    originally 
estimated.       For    this   reason,    and  because   of    the   limited 
time   available   for    this    investigation,    data  from    the   runs 
performed   with   a   gated    sine-wave   input   were   not   analyzed. 
The   following   analysis   pertains    to    the    step-function   input 
data   only. 

Recall    that   the    simple   theory   in   Section   2   based   on 
the    trigometric   approximations   for    the   Bessel   functions 
suggested    that   if   values    of    t^y,    are   plotted   versus    ^^  ,    a 
straight   line    should   result.       This   relationship   would   offer 
advantages    in   data   reduction   because   a   family   of   axis 
crossings    could  be   summarized  by   the    slope   of    a    straight 
line.       For    this    reason,    graphs    which   plotted   the   observed 
and    computed   values    of    Lf^   versus    the    computed   values    of 

Mf*y\      were    made   for    each   run.       Figures   4*2   -   4*4   are 
representative    examples    of    these   graphs.       The    solid   line    of 
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uni 


t   slope   represents    the    computed   values    of    (         .       The 


circles    are    the    computed   values    of    *£*"       ,    and   the    crosses    are 
the    observed   values    of      C^    .       It   is    clear    that    each    set   of 
points    can  be  joined  by   a    straight   line    which   quantitatively 
describes   by   its    slope    the   predicted   and   observed   shift   in   time 
of    the   layer    solution   axis    crossings   as    compared   with    those   of 
the   duct. 

A.  graph   was    then    made   which   plotted   the   predicted   and 
observed   values    of    "^r~~     versus    c/a.(o  n    (Fig.    4.5).       The    simple 
theory   would   suggest   plotting    —22   versus   a/r    (c/a<^j    )     ,    how- 
ever,    r    was   fixed   for   all   experiments    and   only   two   values    of 
a    were   used.       Furthermore   it   is    the    validity   of    the   layer 
solutions    for   various    values    of    c/acon   which   is   the    topic   of   in- 
vestigation here,    rather    than    the   validity   of    Eq.    2.53* 

Analysis   of    the   amplitude   peak   data    was    also    made   for 
step   input    solutions.       The    experimental   amplitude   peaks    read 
from    each   waveform    photograph    were   first   normalized    to    the 
value   of    the   highest   peak   in   the   waveform.       The   predicted 
peaks    were   normalized   in   the    course   of    the    computer   program. 
The   ratio    of    these   normalized   quantities    was    then   plotted   as 
a   function   of    the    delayed   time    of    the   predicted   amplitude 
peaks   for    each   run.       Exact   agreement   in   the    shape   of    the 
transient    waveform   between   experiment   and   theory   would 
result   in  a   line   of    zero    slope.       Figure   4«6    shows    plots    for 
three   values    of    c/ac^n. 
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c/au>K=  0.35 
r  =  61.5  cm 


a  =  2.5  cm 
d  =  2.80  cm 
temp  =  20*  0  • 


Experimental  and  Computed  7^    vs.  £*,  Particle  Displace- 
ment tfaveform  Generated  by  a  Step  Input  in  Velocity. 

Figure  4.2 
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c/aun  =  1.00 

r  -  61.5  cm 


a  =  1.0  cm 
d  =  3 • 14  cm 
temp  =  20  G. 

Experimental  and  Computed  ^  vs.^  Particle  Displace- 
ment waveform  Generated  by  a  Step  Input  in  Velocity. 

Figure  4.3 
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c/au>n  =  3.00 
r  =  61.5  cm 


a  =  1.0  cm 
d  =  9.42  cm 
temp  =  20*  C. 


Experimental  and  Computed  7^,  vs.?^  Particle  Displace- 
ment Waveform  Generated  by  a  Step  Input  in  Velocity. 

Figure  4.4 
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Error   Analysis 

An   analysis    of    the   assumption   that   the    semi-cylindrical 
source   and    steel  backplate    combination    could  be   used   in   lieu 

of   a    cylindrical    source   was    necessary.       The    complex 

( 19) 

impedance    of    the   plate    "was    evaluated  from   relations    in  for 

the   range    of   frequencies    contained   in   the   received   wave- 

J 

forms.       The   impedance    mis-match    was    sufficiently   large   to 
justify   assuming   that   the   plate   was   rigid.       The    directivity 
of    each    transmitter    when    mounted   on   the    steel   plate   was 
investigated   in   order    to    compare   the   resulting   radiation 
pattern   with    the   uniform    radiation   assumed   in   the    theoretical 
deveopment.       Directivity   patterns    are    shown   in   Figs.    i+.7 
and   If.  8.       In   order    to    determine   the   effect   of    the   deviation 
of    the   actual   patterns   from    the   theoretical   patterns t    each 
transmitter    was    clamped   with   foam   rubber    so    that    only   a 
90    degree   arc    in   the    center   of    the    transmitter    was    radiating 
sound.       Except   for    small    change  s    in   amplitude,    no 
differences    were   observed   in   the   received   waveforms    as 
compared   to    those   received   for   radiation   over    the   full    180 
degrees    of    the    transmitter.       It   was    therefore   assumed 
that   radiation   from    the    transmitters    could  be   assumed 
uniform   and    that   the    transmitter-plate    combination   could 
be   used   to    approximate   a    cylindrical   source. 
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The   rise    rime   of    the   experimental   system    was  determined 

to   be   abour.   0.5  jisec.       Therefore,    instead   of   being  a   perfect 

step-function  input,    the   velocity   at   the    source    was  initially 
described  by 

\4a>  -    1  -  * 

where  <*   -   2    x    10      see"         is   the   reciprocal   of    the   rise 

time. 

The   group   velocity   of    the   waveform    is    given  by 

Under    the    condition  Cu>^  <^^    ,    and   with   ~  <*-    ~|"\         Eq.    4,1 

can  be   written   as 

:£.   ^     /-  i  (<"/")* 

By  rearranging   Eq.    l+.l9    the   value   of    delayed   time   at   which 
a   frequency     uj,     will   appear   in   the   received   waveform    can  be 
expressed   as 

It   would  be    expected   that   the   amplitudes    of   frequencies 
greater   than      o^      contained   in   the   initial   portion   of    the 
received   waveforms   would  be   depleted   up   to   a   time    given 

t'  -   ?teO^  u.2> 

and   that   amplitude    distortion   of    the    waveforms    would  be 
significant   up   to    this    time.       However,    substitution   of    the 
observed      o(       and   the   largest   value   of  ui  n   experimentally 
encountered   into    Eq„    4«2,    shows    that    distortion   would  be 

64 


significant   at    the    most   for   only    1.5/JLsec*       Therefore    the 
rise   time   of    the    system    should  have   no    great   effect   on 
the   received   "waveforms. 

It   was    necessary   to   investigate    the    characteristics    of 
the    transmitters    to    see   if    the   velocity   of    the   face 
remained   proportional   to    the   step   input   voltage   over   the 
interval   of    time    that    the   received   waveforms    -were    studied. 
For   times    much   greater   than   the   rise   time,    the   mass   of   the 
diaphram    can  be   ignored   provided   its   velocity   remains 
reasonably   constant   over    the   time   interval   of   interest.      The 

resulting   equivalent    circuit   of    the   transmitter   with   a   source 

(  17) 
resistance   added  is    shown   in   Fig.    k-9-      Writing   loop 

equations   for    the    equivalent    circuit    gives 

«  "  *-  Rs  -  £  i^*  +£fa  V*   Jt  =  ° 

"which    enables   us    to    "write   the   Laplace   transform    of    the 
diaphram   velocity   as 


O     is    determined   from    the    electro-mechanical   relationships 

(  17) 
of    the   transmitter   and   is    defined   as         ' 
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Nominal   equivalent    circuit   values   for    the    1    cm    radius 
transmitter    with   a   20    volt    step   input   applied   are    shown 
in   Fig8    4'9.       Substituting   these   values    into    Eq.    4«3 
and   taking   the   inverse   transform   yields 

-4       -3&>D"£ 


t  ^    3.4  x 


Since   the    mass   of    the    transmitter    diaphram   is   ignored   in 
this    approximation,    the   particle    velocity   is   a    maximum 
immediately   after    the   application   of    the    step   input   and   then 
decays    exponentially.       Substituting  t-j60  /Asec    (the   longest   in- 
terval of  delayed    time    over    which   received   waveforms    were 
observed)    into    Eq.       4»3      shows    that   the   velocity   has 
decayed   to    88        percent   of    its    maximum    value   at   the    end 
of    this    time.       The    tailing   off    of    the   amplitude   peaks    of 
the   received   waveforms,    as    shown   in  Fig.    4.6,    is    consistent 
with    the   decay   in   velocity   of    the   transmitter    diaphram. 

Uncertainties    in   a,    c,    r ?  CO     9    and  7?^    were    the 
possible    sources    of    error   in   this   investigation.       The    effective 
radius    of    each    transmitter    was    measured   to    the   nearest 
1/2    mm,    which    could  have    introduced   an  uncertainty    of    +    2 
percent   in   the   value   of    a   for   the   2.5    cm   radius    transmitter, 
and   +    5   percent   in   the   value   of    a   for    the    1.0    cm    radius 
transmitter.       Measurement   of    the    temperature   of    water 
to    the   nearest   degree   restricted   the    (systematic)    error   in 
c    to    within   +0.5    percent.       The    method   of    determining   T, 
the    time-of -flight   of    the    signal   from    transmitter    to 
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Equivalent  circuit  of  transmitter 
Figure  4.9 
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receiver,    gave   a   possible    systematic    error    of   +    0.5   percent, 
which   in   turn   resulted   in   an   estimated   error   in   r    of    +    0.5 
percent. 

Two    sources    of    error    were   possible   in   the    determina- 
tion  of(j        by   excitation   of    the   tank   as    a   pressure-release 
cavity.       Observed   irregularities    of    +    1/2+    mm    in   the  bottom 
surface    of    the    tank    could   introduce   an   error   in   the   local 
effective   depth    of    the   layer.       Since    experiments    were 
always    conducted   in   the   same    portion   of    the   tank,    a 

systematic    error   in  UJ       would   result.       This    error    could 

n 

range   from    +    0.3    percent'  (for   the   9.2+2    cm    depth)    to   +    1 

percent    (for    the    2.80    cm   depth).       The   random    error 

possible   in   determining   the   frequency   of    the   resonance   peaks 

of    the    standing   wave    spectrum   proved   to  be   insignificant 

relative    to    the   possible    systematic    error. 

The    experimental   values    of  T'T     could  be   influenced  by 

**\ 

systematic    errors    on   the    order   of    +    1   percent    caused  by 

-  "i. 
nonlinearity   and   miscalibration   of    the   oscilloscope   time-base. 

As    mentioned   previously   in   this    section,    axis    crossing   times 
were   read  from    a    series    of    photographs    with   different 
time    scales.       This    method   was   used   so    that   the   initial 
portion   of    the   received   signal   would   always   be   at   the   left- 
hand    edge    of    the   picture   and    the   axis    crossings    of    interest 
would  generally   occur    in   the   right   half   plane   of    the   display. 
This   procedure   resulted   in   a   random    error   of   +   0.5    percent 
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on   all   time    scales    since   the   display   could  be   read   to   within 
l/if    scale   division   out   of    50    divisions. 

The    errors    in   the    computed   values    of  £  .^   "were   influenced 
predominately   by   the    error   in  CO    n,    since      CO   n  appears 
in   the   argument   of    the   Bessel   functions    (    CO^y"*-  —  T*"     )      from 
which    the    C ry\   's       are    computed.       The   values    of       6^        were 
determined   when   the   argument    equaled   any   one    of    a   set   of 
constants      £\       which   resulted   in   a    zero   value   for    the 
amplitude   of    the   particle-displacement   solution.      Thus    the 
argument    can  be   written   as 


and   if    we   have     ^"      <<„     |  ,    th 


*v> 


en 


«ouV5?  r*  *  c 


and 

^    -  *■<"/*?  XT    • 

Therefore   the   +    1   percent   error   in  (0        could  lead   to   +   2 

—        *  n  — 

percent   error   in   the    computed  values    of    ^T^ 

The    errors    in   the    slopes    of    the    curves   fitted   to    the 
experimental    t's'm   Figs.    i+.2    -   2+..J+   were   determined    mainly 
by   the    systematic    errors   in   these   values   of     2"^i     because 
the   random    errors    were    minimized   by  fitting   a    straight   line 
through    the   points.       Therefore    the    error   in   each    slope 
was    estimated    to   be   +    1    percent.     There   was   a    systematic 
error   of   +   2   percent   in   the    computed   values   of     '£"'        ;    hence 
the    error    in    the    slope   of    the    computed      ?~         versus      f 
lines    was   +   2    percent. 
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The   errors   in    c,    a,    and  to      being   independent    can  be 
combined  by   the   root-mean-square   law    to   yield   the    error 

in   c/a  LP  „ .       The   percent   error    in    c/a  u>       varied   according 

'  n  F  n 

to    the   values    of    a   and  &Jn   involved.       Table   i+.l    shows    the 
result   of    this   combination   for   representative   values    of 
c/a  LP     .       Error   flags    on   Figs.    i+.2   -   i+.5   indicate    the 
relative    magnitudes    of    errors    in   the    experimental   and 
computed   values    of     <2^/K1  and   slopes,    and   the   relative 
magnitude   of    error   in   c/a  u>      . 


Table   i+.l 


Tabulated    Errors   for    Quantities    Displayed   on   Figures    i+ .  2—2+ .  5 


Percent    E 

rror 

Exp. 

Comp. 

c/a<on 

c 

a 

a)n 

c/acon 

slope 

slope 

0.35 

+0.5 

+2.0 

+1.0 

+2.3 

+  1.0 

+2.0 

1.00 

0.5 

2.0 

0.3 

2.1 

1.0 

2.0 

1.00 

0.5 

5.0 

0.8 

5.1 

1.0 

2.0 

3.00 

0.5 

5.0 

0.3 

5.0 

1.0 

2.0 
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5.         RESULTS    AND    CONCLUSIONS 

The   investigation  has    shown   that   for    c/a  ^n<^l> 
the   predicted   and    experimentally   observed   shift   in   axis 
crossings    agreed   to    within  +   2   percent.      Where    the 
restriction    c/a  u>        <   1    was   violated,    the   theory  failed   to 
predict    the   observed   shift   in   axis    crossings.      As    c/a  ui 
was    increased    the   failure   of    the   theory  became    more 
pronounced. 

The    correlation  between   the   predicted   and   observed 
amplitudes   was   less   favorable.       The   divergence  between 
the   predicted   and   observed   amplitudes   in   the   initial  portion 
of    the   received    signal   was   unexplained.       The   decrease   in 
the   observed   amplitudes   in   the   latter   portion   of    each 
received   signal   is   not   inconsistent   with    the   roll-off    of    the 
step-function   input   in   velocity. 

Therefore,    with   the    exception   of   amplitudes,    we    can 
say   that    the    theoretical    solutions    correctly   described   the 
propagation   of    simple    transients   in   an   isovelocity   layer    with 
perfectly-reflecting   boundaries   as   long  as    the   parameter 
c/a  do        was    —  1;    and   that,    although    the   theory  failed   for 
values    of    c/a  ui       >    1,    it   did    correctly   predict   the    trend 
that   the   observed    waveforms    followed. 
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APPENDIX   A 
A-l :     Summary   of    the   Propagation   of    Transients    in  Waveguides 

Results   for    the   propagation   of    simple   transients    (step- 
function   and   gated   sine-wave   input)    in   ducts    which  have 
appeared   in   the   literature         '  and   which   will  be    of    value 

in   this   research   are   presented  below    along   with    some   addi- 
tional  Laplace   transform    pairs    that   are   needed. 


1 


Figure  A.l 
Assume   a    waveguide   of    uniform    cross-section,    infinite 
length,    and   perfectly-reflecting  boundaries,    which   is   filled 
with   a   lossless   fluid   and   excited  by   a   velocity    souce   at    z    =   0 
See   Fig.    A.l      The   wave    equation  for    the   acoustic   velocity 
potential     (±)       can  be   written   in    cylindrical    coordinates   as 
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where   r   is    the   radial   direction,    z   is    the   longitudinal   direction* 
and   radially   symmetric   processes   have   been   assumed. 
The    separation   of    variables 


leads    to    the    solution   R      =   *3    (  Co   r/c)    where   «J_    represents 

n  o        n  ° 

the    zeroeth   order   Bessel  function   of    the   first  kind,    and 
the   partial   differential   equation 


jf-^-z^^lt)^  (A'3) 


which   describes    the   propagation   in   the      z         direction. 
The   value   of    the    separation    constant    ^       depends   on  the 
nature   of    the   boundary  at   r   =   a    (the   wall   of    the   duct)    and 
the    cross-sectional   area   of    the   tube.      Because   the   transverse 
properties    of    the   duct   affect   only  60       and   not   the  form    of 
the    solution   of    Eq.    A.  3,    all   functions    of   r   have   been 
suppressed.       They   can  be   reintroduced   with   no   difficulty 
when   desired. 

The   Laplace   transform    of    Eq.    A.  3    -with   quiescent   initial 
conditions    is 

1^-        sVoC"   g  (A. h) 

V 
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where    the  bar    over    the   function  indicates    the   Laplace 
transform    of    that   function,    and    0,        is   a  fragment   of    the 
velocity   potential    containing   z    and   t.       Solving    Eq.    A.  4   with 
the  boundary    condition     ^      =  ^    (0,t)    gives 

0„    ~    \    expC-TvTi^^]  <A-5) 

"where   T    =    z/c   is    the   time   of   flight   of    the   signal   from    the 
source    to    the   receiver. 

The   relation     U.~n   =  Rn    zj^n  is   used   to   obtain   the 

Laplace   transform    of    the   fragment   of    the   longitudinal 
particle   velocity: 


c 

Application   of    the   boundary   condition  on  velocity   at   the 

source   v    (0,t)    =  V    (t)    casts    Eq.    A.6    into    the   following  form 


Of      ^V0   expj^T^^]  (A. 7) 


The   relation 


0o    =     -CX  {A'8) 
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is    evident  from    examination   of    Eqs .    A.  6    and  A.  7.       If 
this   value   of   0      is    substituted   into    Eq.    A.  5,    the   Laplace 
transform   of   0n   in   terms   of    the   velocity  input   becomes 

\  --  - <&  g*pEWCT7_  (A<9) 

The   fragment   of   the   longitudinal   acoustic   pressure   is 
obtained  from    the   velocity   potential  by   means   of    the 
expression    j^  n  =  ^>/p  /r)  =  —0&0r>      .      This    corresponds    to 
p        =   -      &  s0      in   the   Laplace   transform   domain.      We 
therefore   have    the   result 

-£   =   fcs7a_     expfWT^vJ  (A-10) 

The  fragment   of   longitudinal   particle   displacement   may 
be    obtained    either   by   integrating   the   fragment   of   longitudinal 
particle   velocity   in   the    time    domain   or   by   multiplying  by    l/s 
in   the   Laplace   transform    domain.      With    the   latter   alternative, 
the   particle   displacement   fragment      \-»t*         can  ^e   written   as 


-    V£  e*p£Tif?R5fl  (a. ii) 
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Equations  A.  7    and   A.  9   -  A.  11   are    the   Laplace   transforms    of 
the   acoustic   variables    of    interest.       The   inverse   transforms 

for   a   step-function   input    (V      =    l/s)    and  a  gated   sine-wave 

—  °° 

input    (V      =    "TT7!t.         )    are   as   follows: 
o  S  "*o 

For   a    step-function   input 


^30  B4eJ*(«R*} 


Where  ARG   -     O)       y       t2-T2 


For   a   gated   sine-wave   input: 
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(A. 12) 


v>n  --  J^Bfi)  v  2  £  (tlf  T  (Ate)  (a. i3) 


(A. 14) 


(A. 15) 


^  ^[U>fc"V«*)  - U. (n^ «M)|  < A« l6> 


a^k%—    U,(n*rtJfte^  +  O^/ur.jAec)  (a.  17) 


^n"  ^rLu.C-ur^ea)- 0,(^*66)1  (A,18) 


(A. 19) 


Where  ARG   =     Ll>  n    V/       t2-T2 

The   LPs   represent   Lommel  functions   of   two 
variables 


A-2:     Development   of  Additional   Laplace   Transform   Pairs 
for   Step   Input   in  Velocity   in   an  Isovelocity   Layer. 

The   transform   pairs    developed   in   this    section   are 
those   needed   for    the   layer    solution  in  addition   to    those 
summarized   in   section  A-l.       F   denotes    the   transform   under 
consideration;    f(t)    is    the   inverse   transform    of    F.      All  f(t)'s 
are    defined   as    zero   for    t       j£      T. 
1. 

?  ^    sV6   CXpfcTl/IViV1]  (A. 20) 
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With  V      =    1/s,    Eq.    A.  20   becomes 


p      -        <fxP  f-TV*ZHA^l  (A.21) 

Equation  A..  21   is    of   the    same   form    as    Eq.    A.  7   with 

_  CJ 

V      =  _  and     U>  =■  C^  .      Therefore   the   inverse 

transform   of    Eq.    A.  7    may  be   used   with    CO » u±n      •       The 
inverse   transform   of    Eq.    A.  7   is 

O^fl^wr*)  ■|-0l('uJ-,  *e6>)  (A. 22) 


where  ARG   =     (A    U  t  -T 

n 


P^ 


and  /UT±  „    W^^V/- 9/A 

The   quantity     c  «   /c   can  be   written  as      y  |    — .  £^w/ui\ 
With       to  =■  t^^  ,    we   have    /W^.  =.  /Wi   -   u>n(t-T). 

This   in   turn   means    that  both   terms   of    Eq.    A. 22   are 

/  -i  i  \ 
equal.       Using   the   relation  , 


the   desired   inversion   is    obtained.      With  \)       =1   and 
ARG   =  k>„     \|t2-T2,    Eq.    A. 21    may  be   rewritten   as 


80 


2. 


(A. 25) 


With   Vq    =    1/s,    Eq.    A. 25   becomes 


"p   _       e*f>[-T\|sViA>] 


S^VsVod^' 


The    identity 


06 


„0\>-H 


(A. 26) 


(A. 27) 


is   used   to   generate   the   desired   inverse   transform.      Squaring 
Eq.    A.  27    gives 


(it  -  h  I  w  feM  • 


(A. 28) 


Letting   x   =   <^n/s   in    Eq.    A. 28   and   rearranging  yields 


<*>  *  *  l<*b&%  ■ 


(A. 29) 


Solving   Eq.    A.  29   for    l/s   and   substituting   the   value   obtained 
into    Eq.    A.  26    results    in 


?.± 


U> 


(A. 30) 
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(8) 
The   inverse   transform    of    Eq.    A.30   is   available        .       There- 
fore  with       ^   +    1   =  ff)     * 

do 

no 

(A. 31) 


*"i  MS)  J» 


2*1 


A-3 :     Inverse   Laplace   Transforms   for   a   Gated   Sine-Wave 
Input   in  Velocity  in  an   Isovelocity   Layer, 

The  following  acoustic   quantities    could   not  be 
completely   solved.      Although   the   expressions   listed  are 
incomplete   solutions,    it   is   worthy  of   note   that   the   leading 
terms   are   the    corresponding   terms    of    the   duct   solutions 
in   Section  A-l.      In  the   equations   listed  J      j      £      denotes 
the   inverse   Laplace   transform   of   the   quantity   in  brackets. 


(A. 32) 


*  f^  ev[ui^^  -  u«  ^v*^  «  &»(if$y&4r^Jl 


(A. 33) 
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where  ARG   =  (^)     lljl.-ri  ' 

Uyt/W^ARfc)  =    Lommel   Function   of    two   variables^    ^»15^. 

■ft  /       iVtfM 
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DISPLACEMENT 
AXIS  CROSSINGS 
TPRI  IN  USEC 


APPENDIX  B 
EXPERIMENTAL  DATA 

DISPLACEMENT 

AMPLITUDES 

(NORMALIZED) 


VELOCITY 

AXIS  CROSSINGS 

TPRI  IN  USEC 


STEP    INPUT    RUN    S3 

C/AWN=0.35  R=61.5CM      A=2.54CM      D=2.80CM      TMP=20DEGC 


0,0 

1.10 

2.85 

5.50 

9.00 

13.3 

18.5 

24.3 

30.8 

38.0 

45.9 

54.2 

63.8 

73.5 

84.0 

96.0 


0.175 
■0.379 

0.442 
■0.592 

0.660 
■0.718 

0.864 
-0.922 

1.000 
-0.898 

0.893 
-0.718 

0.718 
■0.592 

0.563 


0.95 
2.35 
4.35 
7.65 
11.3 
16.5 
21.7 
28.2 
34.3 
42.4 
50.3 
59.2 
68.2 
79.5 
91.0 


STEP  INPUT  RUN  S12 

C/AWN=0.52     R=61.5CM   A=2.54CM   D=4.15CM   TMP=20DEGC 


0.0 

2.50 

5.10 

10.7 

17.7 

26.1 

36.6 

48.4 

61.3 

76.0 

91.7 

108.5 

126.0 

145.0 

164.0 

184.5 

0.198 
-0.396 

0.636 
-0.833 

1.000 
-0.972 

0.905 

-0.860 
0.838 

-0.838 
0.721 

-0.694 
0.622 

-0.613 
0.532 

-0.477 


NOT  RECORDED 
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DISPLACEMENT         DISPLACEMENT         VELOCITY 

AXIS  CROSSINGS       AMPLITUDES  AXIS  CROSSINGS 

TPRI  IN  USEC         (NORMALIZED)         TPRI  IN  USEC 


STEP  INPUT  RUN  S6 

C/AWN=0.69     R=61.5CM   A=1.00CM   D=5.50CM   TMP=20DEGC 

0.0  0,414  2.20 

3.50  -0.687  6.70 

9.60  0.990  14.4 

18.4  -1.000  23.2 

29.8  0.919  38.6 

45.8  -0.909  52.5 

64.0  0.828  73.5 

85.0  .  -0.854  93.0 

106.0  0.798  119.0 

131.0  -0.727  140.0 

168.0  0.768  169.0 

183.0  -0.687  193.0 

212.0  0.687 

240.0  -0.657 
268.0  0.667 

298.0  -0.581 
329.0  0.616 

360.0 


STEP  INPUT  RUN  S13 

C/AWN=0.85     R=61.5CM   A=2.54CM   D=6.78CM   TMP=21DEGC 

0.0  0.357                NOT  RECORDED 

4.20  -0.697 

13.6  0.871 

26.3  -0.874 

44.5  0.994 

65.9  -0.931 

90.2  0.777 

118.0  -0.834 

147.5  1.000 

179.5  -0.943 

212.5  0.783 
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DISPLACEMENT         DISPLACEMENT         VELOCITY 

AXIS  CROSSINGS       AMPLITUDES  AXIS  CROSSINGS 

TPRI  IN  USEC         (NORMALIZED)         TPRI  IN  USEC 


STEP  INPUT  RUN  S14 

C/AWN=1.00     R=61.5CM   A=2.54CM   D=7.97CM   TMP=21DEGC 

0,0  0.538                NOT  RECORDED 

7.00  -0.813 

19.0  0.869 

36.0  -0.875 

59.5  0.943 

89.2  -0.956 

120.0  1.000 

156.0  -0.863 

190.0  0.769 

235.0  -0.619 


STEP  INPUT  RUN  S4 

C/AWN=1.00     R*61.5CM   A=1.00CM   D=3.14CM   TMP=20DEGC 

0.0  0.176  0.80 

1.25  -0.412  2.44 

3.45  0.559  4.90 

6.40  -0.598  8.75 

10.5  0.696  13.3 

16.0  -0.814  28.9 

21.9  0.961  25.2 

28.2  -1.000  32.4 

35.7  1.000  40.5 

45.2  -0.971  49.2 

53.5  0.922  58.0 

63.5  -0.863  68.5 

74.0  0.784  79.0 

85.0  -0.706  91.5 

98.5  0.667 
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DISPLACEMENT         DISPLACEMENT         VELOCITY 

AXIS  CROSSINGS       AMPLITUDES  AXIS  CROSSINGS 

TPRI  IN  USEC         (NORMALIZED)  TPRI  IN  USEC 


STEP  INPUT  RUN  S9 

C/AWN=1.50     R=61.5CM   A=1.00CM   D=4.71CM   TMP=21DEGC 

0.0  0.307  1.30 

2.30  -0.586  4.20 

6.50  0.793  9.20 

12.8  -0.923  16.5 

20.6  1.000  25.0 

30.6  -0.964  35.7 

42.5  0.979  49.0 

57.0  -0.943  65.0 

71.5  0.821  79.0 

88.0  -0.686  98.0 

107.0  0.750  116.0 

125.0  -0.700  135.0 

145.0  0.593  154.0 

169.0  -0.543  178.0 

191.0  0.521  202.0 

217.0  -0.543  225.0 

239.0  0.564  246.0 

262.0  -0.479  269.0 

286.0  0.486  297.0 
312.0 


STEP  INPUT  RUN  S5 

C/AWN=2.00     R=61.5CM   A=1.00CM   D=6.28CM   TMP=20DEGC 

0.0  0.351  2.10 

3.80  -0.649  7.60 

11.3  1.000  15.5 
21.1  -1.000                                        26.2 

34.4  0.979  41.0 

52.5  -0.918  59.5 
72.5  0.907                                         79.5 
95.5  -0.701  101.0 
119.0  0.731  127.0 
145.0  -0.639  155.0 
176.0  0.639 

203.0  -0.639 

235.0  0.639 

265.0  -0.598 

295.0  0.495 
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DISPLACEMENT         DISPLACEMENT         VELOCITY 

AXIS  CROSSINGS       AMPLITUDES  AXIS  CROSSINGS 

TPRI  IN  USEC         (NORMALIZED)         TPRI  IN  USEC 


STEP  INPUT  RUN  S7 
C/AWN*3.00     R=61.5CM 

0,0 

7.50 

21.0 

40.5 

67.0 

100.0 

132.0 

184.0 

222.0 

272.0 

312.0 


A-1.00CM   D=9. 

42CM   TMP=20DEGC 

0.432 

4.50 

0.892 

15.0 

1.000 

31.0 

0.826 

53.0 

0.822 

84.5 

0.784 

115.0 

0.526 

150.0 

0.427 

210.0 

0.376 

0.319 

0.639 
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VELOCITY  VELOCITY 

AXIS  CROSSINGS       AMPLITUDES 
TPRI  IN  USEC         (NORMALIZED) 


GATED  SINE  INPUT  RUN  Gl  W=583.3KHZ 

C/AWN»1.00     R=6l.5CM   A=2.54CM   D=7.97CM   TMP=20DEG 

3.80  0.279 

9.60  -0.669 

15.9  1.000 

21.5  -0.961 

26.3  0.747 

31.9  -0.870 

37.0  0.851 

42.2  -0.792 

47.9  0.909 

53.0  -0.870 

58.5  0.857 

63.5  -0.929 

69.0  0.864 

75.0  -0.909 

79.5  0.870 

84.5  -0.877 

89.5  0.922 

95.0  -0.864 

100.5  0.916 

105.0  -0.877 

115.5  0.890 

120.5  -0.883 

125.5  0.857 

131.0  -0.896 

136.5  0.864 

141.5  -0.877 

147.0  0.883 

152.5  -0.857 

157.5  0.896 

164.0  -0.831 

172.5 
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